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Abstract. This article presents a phenomenological dynamic phase transition 
theory for ferromagnetism, leading to a precise description of the dynamic tran- 
sitions, and to a physical predication on the spontaneous magnetization. The 
analysis also suggests asymmetry of fluctuations in both the ferromagnetism 
and the PVT systems. 



1. Introduction 

Classically, phase transitions are classified by the Ehrenfest classification scheme, 
based on the lowest derivative of the free energy that is discontinuous at the transi- 
tion. For ferromagnetic systems, it has been observed that the magnetization, which 
is the first derivative of the free energy with the applied magnetic field strength, 
increases continuously from zero as the temperature is lowered below the Curie tem- 
perature, and the magnetic susceptibility, the second derivative of the free energy 
with the field, changes discontinuously. Hence, the ferromagnetic phase transition 
in materials such as iron is regarded as a second order phase transition. However, 
a theoretical understanding of the transition is still lacking. The main objective 
of this article is to provide theoretical approach to dynamic phase transitions for 
ferromagnetic systems. 

For the classical GL free energy, although both the steady state and time- 
dependent models provide some results in agreement with experiments, there are 
obvious discrepancies on both susceptibility and spontaneous magnetization. Hence 
a revised GL free energy is proposed and analyzed, leading to a precise description 
of the dynamic transitions, and to a physical predication on the spontaneous mag- 
netization. 

The analysis is based on the recently developed dynamic transition theory by 
the authors, together with a new dynamic classification scheme, which classifies 
phase transitions into three categories: Type-I, Type-II and Type-Ill, correspond- 
ing mathematically to continuous, jump, and mixed transitions, respectively; see 
Section 2 and the Appendix as well as two recent books by the authors [U [S] for 
details. 

We remark also that the analysis leads naturally to a physical conjecture on 
asymmetry of fiuctuations, which appears in both the ferromagnetic system studied 
in this article and in PVT systems studied in [2]. 
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This article is organized as follows. In Section 2, we review the dynamic classifi- 
cation scheme and the new time-dependent Ginzburg-Landau model for equilibrium 
phase transitions. Section 3 deals with the dynamic transition based on the classi- 
cal Ginzburg-Landau energy, and Section 4 addresses the dynamic transition theory 
using a revised Ginzburg-Landau energy. Physical conclusions are given in Section 
5, and dynamic transition theory is recapitulated in the Appendix. 



2. General Principles of Phase Transition Dynamics 

In this section, we recapitulate the new dynamic phase transition classification 
scheme to classify phase transitions into three categories: Type-I, Type-II and 
Type-Ill, corresponding mathematically to continuous, jump, and mixed transi- 
tions, respectively. 

Then we recall a new time-dependent Ginzburg-Landau theory for modeling 
equilibrium phase transitions in statistical physics, derived based on the le Chatelier 
principle and some mathematical insights on pseudo-gradient systems. 

Both the classification scheme and the Ginzburg-Landau theory was developed 
recently by the authors, and we refer interested readers to [TJ|3J|2] for details. 

2.1. Dynamic classification scheme. In sciences, nonlinear dissipative systems 
are generally governed by differential equations, which can be expressed in the 
following abstract form 

(2.1) -^^=L,u + G{u,X), 

u{0) = (fi, 

where u : [0, oo) — > X is the unknown function, A e {N > 1) is the control 
parameter, X and Xi are two Banach spaces with Xi C X being a dense and 
compact inclusion, L\ = —A + B\ and G'(-, X) : Xi ^ X are C"'(r > 0) mappings 
depending continuously on A, L\ : Xi X \s a, sectorial operator, and 

A : Xi ^ X a linear homeomorphism, 

(2.2) 

B\ : Xi ^ X a, linear compact operator. 

In following, we introduce some basic and universal concepts in nonlinear sci- 
ences. 



First, a state of the system (2.11 at A is usually referred to as a compact invariant 



set Yi\. In many applications, Yj\ is a singular point or a periodic orbit. A state 



Ti\ of (2.11 is stable if Ti\ is an attractor, otherwise Y^x is called unstable. 



Second, we say that the system (2.1 1 has a phase transition from a state Sa at 
A = Ao if Sa is stable on A < Aq (or on A > Aq) and is unstable on A > Aq (or on 
A < Aq). The critical parameter Aq is called a critical point. In other words, the 
phase transition corresponds to an exchange of stable states. 

The concept of phase transition originates from the statistical physics and ther- 
modynamics. In physics and chemistry, "phase" means the homogeneous part in 
a heterogeneous system. However, here the so called phase means the stable state 
in the systems of nonlinear sciences including physics, chemistry, biology, ecology, 
economics, fluid dynamics and geophysical fluid dynamics, etc. Hence, here the 
content of phase transition has been endowed with more general significance. In 
fact, the phase transition dynamics introduced here can be applied to a wide variety 
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of topics involving the universal critical phenomena of state changes in nature in a 
unified mathematical viewpoint and manner. 



Third, if the system (2.1 1 possesses the gradient-type structure, then the phase 
transitions are called equilibrium phase transition; otherwise they are called the 
non-equilibrium phase transitions. 

Fourth, classically, there are several ways to classify phase transitions. The one 
most used is the Ehrenfest classification scheme, which groups phase transitions 
based on the degree of non-analyticity involved. First order phase transitions are 
also called discontinuous, and higher order phase transitions {n > 1) are called 
continuous. 

Here we introduce the following notion of dynamic classification scheme: 



Definition 2.1. Let Xq E M.^ be a critical point of {2.1), and (2.1) undergo a 
transition from state Yi\ to Yi\. There are three types of phase transitions for (2.1) 
at X = Xq, depending on their dynamic properties: continuous, jump, and mixed as 
given in Theorem \A.l\ which are called Type-I, Type- II and Type-Ill respectively. 

The main characteristics of Type-II phase transitions is that there is a gap be- 
tween Sj^ and at the critical point Aq. In thermodynamics, the metastable states 
correspond in general to the super-heated or super-cooled states, which have been 
found in many physical phenomena. In particular, Type-II phase transitions are 
always accompanied with the latent heat to occur. 

In a Type-I phase transition, two states and meet at Aq, i.e., the system 
undergoes a continuous transition from to 

In a Type-Ill phase transition, there are at least two different stable states 
and Sg at Aq, and system undergoes a continuous transition to E2 or a jump 
transition to S3, depending on the fluctuations. 

It is clear that a Type-II phase transition of gradient-type systems must be 
discontinuous or the zero order because there is a gap between and I]]^ (2 < j < 
K). For a Type-I phase transition, the energy is continuous, and consequently, it 
is an n-th order transition in the Ehrenfest sense for some n > 2. A Type-Ill phase 
transition is indefinite, for the transition from to it may be continuous, i.e., 
dF~ /dX — dF2 /dX, and for the transitions from T,\ to it may be discontinuous: 
dF'/dX ^ dF+/dX at A = Aq. 



2.2. Time-dependent Ginz burg- Landau models for equilibrium phase tran- 
sitions. In this subsection, we introduce the time-dependent Ginzburg-Landau 
model for equilibrium phase transitions. 

We start with thermodynamic potentials and the Ginzburg-Landau free energy. 
As we know, four thermodynamic potentials- internal energy, the enthalpy, the 
Helmholtz free energy and the Gibbs free energy-are useful in the chemical ther- 
modynamics of reactions and non-cyclic processes. 

Consider a thermal system, its order parameter u changes in il C M" (1 < n < 3). 
In this situation, the free energy of this system is of the form 



/" 1 

(2.3) niu,X)=na+ [-^/i.|Vu,|2-f 5(u,Vm,A) 



dx 
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where > 3 is an integer, u ~ [ui, • • • , Um), ^J■i — ^J■iW > 0, and g{u, Vm, A) is a 
C"'(r > 2) function of (w, Vu) with the Taylor expansion 

N 

(2.4) g{u, Vu, A) = ^ a^JkU^DjUk + ^ a/u^ + o(|u|^) - /X, 

|/|=i 

where / = (ii,--- ,im),'ik > are integer, |/| = X)feLi*fc> the coefficients a^fc 
and a I continuously depend on A, which are determined by the concrete physical 
problem, = ■ ■ ■ ul^ and fX the generalized work. 

Thus, the study of thermal equilibrium phase transition for the static situation 
is referred to the steady state bifurcation of the system of elliptic equations 

du 

-T^ldn^O, (orM|an = 0), 

where S/Su is the variational derivative. 

A thermal system is controled by some parameter A. When A is for from the 
critical point Aq the system lies on a stable equilibrium state Si, and when A 
reaches or exceeds Aq the state Si becomes unstable, and meanwhile the system 
will undergo a transition from Si to another stable state S2. The basic principle 
is that there often exists fluctuations in the system leading to a deviation from the 
equilibrium states, and the phase transition process is a dynamical behavior, which 
should be described by a time-dependent equation. 

To derive a general time-dependent model, first we recall that the classical le 
Chatelier principle amounts to saying that for a stable equilibrium state of a system 
S, when the system deviates from S by a small perturbation or fluctuation, there 
will be a resuming force to retore this system to return to the stable state S. Second, 
we know that a stable equilibrium state of a thermal system must be the minimal 
value point of the thermodynamic potential. 

By the mathematical characterization of gradient systems and the le Chatelier 
principle, for a system with thermodynamic potential 7i(u, A), the governing equa- 
tions are essentially determined by the functional 7i(u, A). When the order param- 
eters (ui, • • • ,Um) are nonconserved variables, i.e., the integers 

Ui{x,t)dx — ai{t) ^ constant. 

n 

then the time-dependent equations are given by 
du. 
~di 

(2.5) du 



= -A— H(w, A) + $,(m, Vw, A), 
at ouj 



Q^idn^O (or = 0), 



u{x, 0) = ip{x), 

for any 1 < « < m, where S/6ui are the variational derivative, Pi > and $i satisfy 



(2.6) 



/ y"$i-^H(u,A)da; = 0. 



The condition (2.6) is required by the Le Chatelier principle. In the concrete 



problem, the terms i'i can be determined by physical laws and (2.6 1 
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When the order parameters are the number density and the system has no ma- 
terial exchange with the external, then Uj (1 < i < Tn) are conserved, i.e.. 



(2.7) 



Uj{x, t)dx ~ constant. 



This conservation law requires a continuous equation 



(2 



9u, 



dt ' 

where Jj(u,X) is the flux of component Uj. In addition, Jj satisfy 
(2.9) J,^-k,\/{fi,-Y,f^,), 

where fii is the chemical potential of component ui, 

S 



(2.10) 



-Ti{u, A) — (t>j{u, Vw, A), 



and <j)j{u,\) is a function depending on the other components Ui (i ^ j). When 



m — 1, i.e., the system consists of two components A and B, this term 
Thus, from (2.8l-(2.10l we obtain the dynamical equations as follows 



0. 



(2.11) 



Ion = 0, 



dAu. 



dt 
du 

n an 
u{x, 0) — ip{x), 



H{u, A) — (/)j(M, Vw, A) 



Ian = 0, 



for I < j < m, where (3j > are constants, satisfy 



(2.12) 



H{u,X)dx = 0. 



If the order parameters (ui,-- - ,Uk) are coupled to the conserved variables 
(itfc+i, • • • , Urn), then the dynamical equations are 

-l3.~n{u,X) + ^,{u,Vu,X), 



(2.13) 



dt 

dt 

dui 
dn 
duj_ 

dn 



PA 



Su 



-Ti.{u, A) — (t)j{u, Vw, A) 



Ian = (or Uilon 



\dn 



0, 



dAuj 
dn 



I an 



0), 
= 0, 



u{x, 0) = (p{x). 

for I < i < k and k + I < j < m. 

The model (2.13) gives a general form of the governing equations to thermody- 
namic phase transitions. Hence, the dynamics of equilibrium phase transition in 
statistic physics is based on the new Ginzburg-Landau formulation (2.13). 

Physically, the initial value condition u(0) = ip in (2.13) stands for the fluctuation 
of system or perturbation from the external. Hence, (p is generally small. However, 
we can not exclude the possibility of a bigger noise ip. 
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From conditions (2.6) and (2.121 it follows that a steady state solution uq of 



(2.13 1 satisfies 



(2.14) 



$,(uo, Vwo, A) = yi<i<k, 
A(^j(uo, Vmo, A) = Vfc+l<j<m. 



Hence a stable equilibrium state must reach the minimal value of thermodynamic 
potential. In fact, uq fulfills 

(3i—H{ua,X) -$i(uo,Vuo,A) = 0, 

nil . 



(2.15) 



du 



-^|ao = (oruj|an = 0), 



Ian = 0, 



Ian 



0, 



I CI i ^ : ^ 

on on 

for 1 < i < fc and fc + 1 < j < m. Multiplying <i>i(uo, Vuq, A) and 0j(uo,Vwo,A) 



on the first and the second equations of (2.151 respectively, and integrating them, 
then we infer from (2.6 1 and ( |2.12[ ) that 



/ V$?(uo,Vwo,A)da: = 0, 
Jn ^ 

I V|V0j(Mo,Vuo,A)|2dx = O, 
"'n . 



which imply that (2.141 holds true 



3. Classical Theory of Ferromagnetism 

A ferromagnetic material consists of lattices containing particles with a magnetic 
moment. When no external field is present and the temperature is above some 
critical value, called the Curie temperature, the magnetic moments are oriented at 
random and there is no net magnetization. However, as the temperature is lowered, 
magnetic interaction energy between lattice sites becomes more important than 
the random thermal energy. Below the Curie temperature, the magnetic moments 
become ordered in the space and a spontanous magnetization appears. The phase 
transition from a paramagnetic to a ferromagnetic system takes place at the Curie 
Temperature. 

The phase diagrams for magnetic systems are given in Figures |3.1j|3.3| In Figure 
|3.1| below the Curie temperature, the magnetization occurs spontaneously, and the 
zero magnetic field H = separates the two possible orientations of magnetization. 
Figure |3.2| provides a sketch of the isotherms of magnetic system, and Figure |3.3| 
gives the magnetization as a function of temperature; see also Reichl [S] and Onuki 
[4] for details. 

Based on the classical Ginzburg-Landau theory, for an isotropic system, the 
Helmholtz free energy can be expressed as 



A(M,r) = Ao(T) + ia2(r)|M|2 



1 



a,{T)\M\' 
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Figure 3.1. Below the Curie point the magnetization occurs 
spontaneously; the curve H = separates the two possible ori- 
entations of magnetization. 




Figure 3.2. A sketch of the isotherms for a magnetic system. 




Figure 3.3. A sketch of the spontaneous magnetization of a mag- 
netic system. 



where Aq (T) is a magnetization-independent contribution to the free energy, | M ^ 
M ■ M, and AI — {Mi, M2, M^) is the magnetization of the system. When 
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external field H is present, the Gibbs free energy is given by 
G(M , T) = A(M, T)-H -M 
= Ao(T) ~H -M 



^a2(T,i7)|M|2 + ia4(r,i/)|M|4 



For small H , a2 and 014 can be considered to be independent of iJ, and near the 
Curie point we have 

a2{T) = ao{T~-T,), a^iT) > 0. 

Usually, G{M, H, T) is called the Ginzburg-Landau free energy. To omit the 
higher order terms than it is known that the equilibrium state M of the 

ferromagnetic system satisfies 

(3.1) —G^aAM\'^M + a2M-H = Q. 

oM 



Thus above the Curie point we obtain from (3.1 ) that 

1 



(3.2) 



M ~ —H, 
0L2 
dM 

X ■ 



1 



dH 02 (T) ao(T-r,)' 

where x is the isothermal susceptibility, which is a scalar because the system is 
isotropic. Below the critical point, for H = {), the magnetization M obeys 

' d\M\ 



IMT,~T) 



(3.3) \M\ 
The heat capacity aX T — T^ is 

(3.4) C{T < Tc) - C{T > Tc) 



dT 



ao 



ai{T,-Ty 



QT2 



1 



'-a4\M\ 



' 2^4 ■ 



We infer then from (3.2l-(3.4) the following classical conclusions for an isotropic 
magnetic system: 

(1) When as external magnetic field is present, a nonzero magnetization exists 
above the Curie point Tc, which has the same direction as the applied field 
U. 

(2) Near the critical point the susceptibility x tends to infinite with the rate 
(T— Tc)""'^, i.e., a very small applied field a,t T = Tc can yield a large effect 
on the magnetization. 

(3) In the absence of an external field (i.e., H = 0), below the critical point a 
spontaneous magnetization M appears, which depends continuously on T 
and tends to zero with the rate {T — Tc)^^'^; namely the transition is of the 
second order. ^ 

(4) The heat capacity at T — Tc has a jump with the gap AC — and the 



jump has the shape of a A, as shown in Figure [3^ 



2a4 
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Qualitatively, part of the above conclusions are in agreement with experimental 
results. 

However these conclusions lead to wrong susceptibility and spontaneous magne- 
tization, whose experimental rates are given by x (T — Tc)^^ with r = 1.3, and 
by M oc (T - Tc)^ with /3 = 1/3. 

Free energy G must be a function of the magnetization A/. Hence the errors are 
originated from the fact that the expression of G in the Ginzburg-Landau theory 
is an approximation. It is difficult to derive a precise formula because G is not 
analytic on \M\, even the dilTerentibility of G on |M| is very low. 

If we study the dynamical properties of ferromagnetic systems by using the 
classical Ginzburg-Landau free energy, we shall see a more serious error when an 
external field is present. 

To see this, the dynamic equation of classical theory is given by 

(3.5) ^ = -a2M ~ a4|MpA/ + H. 

For simplicity, we take H = (/i, 0, 0) with /i > 0, it is equivalent that we take the 
xi-axis in the direction of H. Then the equation 

ai\M\^M + aaM -H = 

has a steady state solution Mq = (too,0,0) for T > 0, which is the magnetization 
induced by H. Make the transformation 

M ^M' + Mn. 



Then, the equation (3.51 is rewritten as (drop the primes) 
dMi 



dt 



= -{a2 + 3a4TOo)Mi - 2aim(^M( - aima\M\^ - a4|MpMi, 



(3.6) = -{a2 + aiml)M2 - 2aimaMiM2 - ai\M\'^M2, 



dt 
dMs 



'{a2 + Q4mo)M3 - 2a4moMiAf3 - a4|MpM3. 



dt 

Comparing the two critical parameter curves 

a2 + 3Q;4mQ ^ 2 ^ Ti = Tc — iaimfj ao^ 
a2 + a4ml = ^ T2 = Tc — a4m^/ao, 



we find that T2 > T^. By Theorem A.l (pU3| has the &-st transition at T = 



where a new magnetization M = (Mi, M2, M3), with M2 ^ and M3 ^ 0, appears. 
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This is unrealistic because any magnetization M of this system must have the same 



direction as H — (/i, 0,0); see Figure 3.1 



In fact, when a magnetic field H is applied on an isotropic system the direction 
of is a favorable one for magnetization. However, in (3.5 1 this point is not man- 
ifested. Therefore, to investigate the phase transition dynamics of ferromagnetic 
systems we need to revise the GL free energy. 

4. Dynamic Transitions in Ferromagnetism 

4.1. Revised Ginzburg-Landau free energy. Let the ferromagnetic system be 
isotropic. When a magnetic field H is present, we introduce a second order sym- 
metric tensor 

A{T,H) = {ai.j{T,H)), Uij ^ aji, 1 < «, j < 3, 
such that A{T, 0) = 0, and A{T, H) has eigenvalues 

Ai = Ai(T,i/), A2 = A3 = 0, with Ai(r,i7) > as 7^ 0, 
and H is the eigenvector of A corresponding to Ai: 

(4.1) AH^XiH (Ai > as 7^ 0, Ai = as = 0). 

It is clear that if we take the coordinate system (xi, a;2, 2:3) with xi-axis in the 
i?-direction, then H = (J?i,0,0) and 

(4.2) A = 



Ai 



























Physically, condition (4.1) means that H is & favorable direction of magnetization 
if we add a term 

-MAM^ = - ^ aijM.Mj 

in the free energy. We also need to consider the nonlinear effect acted by H. To 
this end we introduce the term — |MpM • H in the free energy. 

Thus when the applied field H may vary in 17 C M" {n = 2, 3), then the GL free 
energy is in the form 



(4.3) G{M,T,H) = Go + l f M\' + a2\M\^ + ^\M\ 

- J2 a^JM,M, - f{T, H, M)M ■ H 



where Go — G(){T) is independent of M and H, and / is a scalar function of T, H 
and M, defined by 

(4.4) f{T,H,M)^2{l + P\M\^), p = P{T,H)>0. 
For a2 and we assume that 

(4.5) a2^ao{H){T-To{H)), ro(0) = T„ ao > 0, 04 > 0. 



By the standard model (2.51, we derive from (4.3 1 and (4.4 1 the following dy- 
namical equations: 

(4.6) ^ ^fiAM, - a2M, + ^ a„M, - a^lMl^M, 

j=i 

-I- /3|Af pffi -t- 2/3(M • iJ)M, + for 1 < i < 3. 
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The boundary condition is given by 



(4.7) 



dM 
dn 



an 



Obviously, if H ~ 0, (4.6 1 coincide with the classical equations. For simplicity, 
hereafter we always take 

(4.8) H = {h, 0, 0) (/i > is a constant). 



When H is constant, in the study of phase transitions of magnetic systems, (4.6) 
can be replaced by a system of ordinary differential equations as follows: 



(4.9) 



dt 

dM2 

dt 
dMs 

dt 



= (Ai - a2)Mi + f3h\Mf + 2(ihMl - ai\M\^Mi + h, 
= + 2l3hMiM2 - a4\MfM2, 

= -Q2M3 + 2l3hMiM3 - a^lM^Ms. 



Equations (4.9 1 have a steady state solution induced by H: 
M* = (too; Oi 0)j with lim mg = 0. 

We see in Figure [3?2| that the magnetization M* has an upper bound. Namely there 
is an Mq such that 



(4.10) 



\M*\ < \Mo\ 
M* ~* Mo 



yn e M-', r > o, 

if h 00. 



To satisfy (|4.10|) it is necessary to assume that the coefficients ag and as in (4.5) 

-00, a^iT, H) — > +00, as /i — + 00. 



possess the properties 
(4.11) ao{H) 



Both conditions ( |4.5| and (4.11) are physical. 

4.2. Dynamic transitions. In this subsection, to illustrate the main ideas, we 
only consider the case where iJ is a constant on il. Therefore we shall study phase 
transition dynamics of the ferromagnetic systems by using equations (4.9 1 for h > 0. 
Analysis for more general case can be carried out in the same fashion, and will be 
reported elsewhere. 

Take the transformation in (4.9) 

(4.12) 



M = M* + M'. 



Then equations (4.9) are rewritten as (drop the primes) 

= /32M2 - 2a2MiM2 - a^\M\^M2, 
P2M3 - 2a2MiM3 - a4|MpM3, 



(4.13) 



dt 
dM2 

dMs 
dt 
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where 

02 — Ct4TnQ — ph, 

(3i — \i + 6phmo — 3Q!4mo — a2, 
/?2 = 2/3/iTOo — Q!4TOq — a2- 
The critical parameter curves Pi — and /32 = are given by 

/3i = Ti = ro(i7) + — (Ai + ephmo - 3a4ml), 

/32 = = To{H) + — (2/3/imo - aiml). 

ao 

It is clear that Ti > T2 provided 
(4.14) Ai > 2mo(a4TOo - 2/3/i), 



for ft, > 0. 



Therefore, under condition (4.14), the equations (4.131 have a transition at T - 
Ti in the space 

(4.15) £:= {(Afi,0,0)| -00 < Afi < 00}. 

More precisely, we have the following transition theorem. 



Theorem 4.1. Assume the condition {4-14^ <^^d, 02 7^ 0. Then has a Type- 

Ill (mixed) transition at T — Ti , and the transition occurs in the space E. The 
phase diagram is as shown in Figure \4^ Moreover we have the following assertions: 
(1) There are two stable equilibrium states near T = Ti, which are given by 
ifT>Ti, 
1 



m; 



2a4 

1 

2q;4 



-3a2 



%al + AaiPi\ ifT<Ti, 



[3a2 



9al 



(2) If T < Ti, is stable in the region < Mi < 00, and is stable in 
-00 < Ml < 0. 

(3) If T > Ti, ~ Q is stable in —6 < Mi < 00 and M^ is stable in 
—00 < Ml < —b, where 

^ — - [3a2 — y uu,2 



2aA 



9ai + 4a4/3i] > 



forT>Ti. 



Proof. It is clear that if (4.141 holds, then 

r <0 ifr>Ti, 

Pi{T) \ =0 ifT-Ti, 
[ >0 ifr<Ti, 

P2{Ti) < 0. 



Hence, by Theorem A.l 



the system (4.131 has a transition at T = Ti. Obviously, 



the space E defined by (|4.15| is the center manifold of (|4.13|) near T = Ti. Hence, 
the rec 

(4.16) 



the reduced equation of (|4.13| on E is expressed as 

PiMi - 3a2M^ - a^M^ 



dMi 
dt 
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Figure 4.1. 



As a2 7^ 0, by Theorem A. 2 we infer from (4.16) that this transition is of type-Ill, 



0. 



and the transition solutions satisfy 

a^M^ + 3a2Mi - pi 

By a direct compute one obtains Assertions (1) and (2). 
The proof is complete. 

5. Physical Conclusions and Remarks 
5.1. Physical predications based on Theorem 



□ 



states of (4.9 1 near T 



4.1 



By ( 4. 12 1 , the stable steady 



Ti are 

M+ = (too + M+,0,0), 

M- = (mo + Mf ,0,0), 
From the physical point of view, it should be 



(5.1) 



M+ > 0, 



Mf < 0, mo + A/j" > 0. 



The condition (5.1 1 requires that 



< 3a2 < aAfriQ 

which is equivalent to 
(5.2) 

where toq > is a solution of the equation 
(5.3) 047710 — SPhrriQ + {a2 — Xi)mo 



f3h < 047710 < 2^^^^ 



{h > 0) 



0, 



Thus, the stable steady states M+ and M~ of (4.9 1 near T = Ti are physical 



near T — T-i 



provided that the cvoeficients a2{T,h),a4{T,h),P{T,h) and Xi(T,h) satisfy (5.2) 



and (5.3 1. In this case the temperature Ti is greater than the Curie temperature 
Tc-. 

Ti{H) > ^ Ti{0) fori? 7^0. 
The two states and M" are mathematically equal, therefore only by Theo- 



rem |4.1| we can not determine the magnetization behaviors of ferromagnetic systems 
near T — Ti. However, we see that the magnetization M"*" is stronger than M~ . 
Physically, it implies that Af+ is favorable in T < Ti, and is in T > Ti. 
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Thus, from Theorem |4.1| there are two possible magnetization behaviors, i.e., two 
magnetization functions: 

^i,{T) = M+{T) = (mo(r) + A/+(r),o,o), 

( M+{T)^{mo{T) + M+iT), 0,0) if T < T,, 
\m-(T) = (mo(T) + Aff(T),0,0) ifT>ri. 

The function fJ.i{T) is continuous on T, as shown in Figure [5TT| a), and its derivative 
is discontinuous at T = Ti : 



,0,0 



3dTa2{Ti)' 

The function 112 {T) has a jump at T = Ti , as shown in Figure |5.1[b) . 



Figure 5.1. (a) The graph of function Hi{T)\ (b) the graph of 
function IJ.2{T). 



On the other hand, by direct computation, the free energies of M+ and M are 
shown to be given by 

(Gimo) ifr>Ti, 

G{M+) = { 1 ^ , 

|G(mo) + -(Af+)2(a2M+-/3i) ifr<Ti, 

G(M-) = G(mo) + ^(Aff )2(a2Mf - 



where mg is the magnetization induced by _ff = (/i, 0,0) satisfying (5.3), and 
/3i(Ti) = 0. It is clear that 

(5.4) G(Af+) > G(Af-) near T = Ti. 

Hence, it follows from ( |5.4[ ) that the magnetization behavior described by ii2{T) 
is prohibited in real world because the free energy can not abruptly increase (or 
decrease) in a temperature decreasing (or increasing) process. Thus, by Theorem 
4.1 and (5.4 1 we can derive the following physical conclusion: 



Physical Conclusion 5.1. When an external field H is present, the magnetization 
Mh(T) of an isotropic ferromagnetic system is continuous on the temperature T , 
and there is a Ti(H) > Tc (T,. the Curie temperature) with Ti{H) —^Tc as H ^ 
such that Mh{T) is not differentiahle atT ~ Ti, whose derivative has a finite jump 

M'H{T^)-M'H{T+)=a>0 (a < 00). 
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Moreover, the graph of Mh{T) — Hi{T) as shown in Figure \5n][ a ), and Mh{T) 
Mq{T) as H ^0 with 



Mo{T) 



ifT>T,, 
where Mg(T) is the spontaneous magnetization (see Figure 3.3). 



5.2. Asymmetry of fluctuations. The above discussions suggest that for ferro- 
magnetic systems, there are two possible phase transition behaviors near a critical 
point, and theoretically each of them has some probability to take place, however 
only one of them can appear in reality. For the ferromagnetic systems we again 
see this situation. This phenomena is also observed in phase transitions for PVT 
systems 

One explanation of such phenomena is that the symmety of fluctuation near a 
critical point is not generally true in equilibrium phase transitions. To make the 
statement more clear, we first introduce some concepts. 

Let G(w, A) be free energy of a thermodynamic system, u — (ui, • • • , u„) be the 
order parameter, and A = (Ai, • • • , Am) the control parameter (n, to > 1). Assume 
that u is defined in the function space L^(ri,IR") and A e M™. Then the space 

X = {(m,A)| u e L2(f],M"),A e M"} 

is called the state space of the system. 

Let (uq, Aq) G a be a stable equilibrium state of the system; namely (uo,Ao) 
is a locally minimal state of G{u, A). We say that the system has a fluctuation at 
(ito, Ao) if it deviates randomly from (uo, Aq) to {u, A) with 

||m-uo|| + |A- Aol > 0. 

In this case, (w. A) is called a state of fluctuation. 

The so called symmetry of fluctuation means that for given r > 0, all states 
(w, A) of fluctuation satisfying 

||u-uo|| + |A-Ao| =r, (u,A)eA, 

have the same probability to appear in real world. Otherwise, we say that the 
fluctuation is asymmetric. 

The observations in both the PVT systems and the ferromagnetic systems strongly 
suggest the following physical conjecture, regarding to the uniqueness of transition 
behaviors. 



Physical Conjecture (Asymmetry of Fluctuations). The symmetry of fluctua- 
tions for general thermodynamic systems may not he universally true. In other 
words, in some systems with multi- equilibrium states, the fluctuations near a criti- 
cal point occur only in one basin of attraction of some equilibrium states, which are 
the ones that can be physically observed. 

Appendix A. Recapitulation of the Dynamic Transition Theory 

In this appendix we recall some basic elements of the dynamic transition theory 
developed by the authors [TJ [3] , which are used to carry out the dynamic transition 
analysis for the ferromagnetism systems in this article. 
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Let X and Xi be two Banach spaces, Xi C X a. compact and dense inclusion. 
In this chapter, we always consider the following nonlinear evolution equations 



(A.l) 



du 
^(0) 



L\u + G(w, A), 



where u : [0,oo) — > X is unknown function, and A S is the system parameter. 

Assume that Lx '■ Xi — > X is a parameterized linear completely continuous field 
depending contiguously on A G M^, which satisfies 



(A.2) 



L\ = —A + B\ a sectorial operator, 

A : Xi ^ X a linear homeomorphism, 

Bx : Xi X a linear compact operator. 



In this case, we can define the fractional order spaces X^^ for a G M^. Then we also 
assume that G{-,X) : Xa ^ X is C^{r > 1) bounded mapping for some < a < 1, 
depending continuously on A G M^, and 



(A.3) 



Giu,X) = oi\\u\\xJ, VAe 



Hereafter we always assume the conditions (A.2 1 and (A.3 1, which represent that 



the system (A.l I has a dissipative structure. 



In the following we introduce the definition of transitions for (A.l ) 



Definition A.l. We say that the system {A.l) has a transition of equilibrium from 
(u, A) = (0, Ao) on A > Aq (or X < Xq) if the following two conditions are satisfied: 

(1) when A < Ao (or X > Xq),u ~ is locally asymptotically stable for {A.l); 
and 

(2) when A > Aq (or X < Xq), there exists a neighborhood U d X of u ~ 
independent of X, such that for any G U \ Tx the solution ux{t,ip) of 
(A.l) satisfies that 



liTasup\\ux{t,ip)\\x > 6{X) > 0, 



lim S{X) > 0, 

where Tx is the stable manifold of u 
A > Ao (or X < Ao). 



0, with codim Tx > ^ in X for 



Obviously, the attractor bifurcation of (A.l) is a type of transition. However, 



bifurcation and transition are two different, but related concepts. Definition |A.1| 



defines the transition of (A.l I from a stable equilibrium point to other states (not 
necessary equilibrium state) . In general, we can define transitions from one attrac- 
tor to another as follows. 



Definition A.2. Let Sa C AT &e an invariant set of (A.l). We say that (A.l) 
has a transition of states from (Sao, Ao) on A > Aq (or X < Xq) if the following 
conditions are satisfied: 

(1) when A < Aq for A > Ao),Sa is a local minimal attractor, and 

(2) when X > Xq (or X < Xq), there exists a neighborhood U d X of 
independent of X such that for any ip Cz U \ (Tx U S^), the solution u(t, ip) 
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of (A.l) satisfies that 

limsupdist{u{t,ip),T,\) > S{X) > 0, 

t — >co 

lim (5(A) = (5 > 0, 

A— »Ao 

where Tx is the stable manifolds ofY,\ with codim > 1. 
Let the eigenvalues (counting multiplicity) of L\ be given by 

{/3,(A)eC I j = l,2,...} 




VI < i < m, 
Vj > 7n + 1 . 

The following theorem is a basic principle of transitions from equilibrium states, 
which provides sufficient conditions and a basic classification for transitions of non- 
linear dissipative systems. This theorem is a direct consequence of the center man- 
ifold theorems and the stable manifold theorems; we omit the proof. 



Theorem A.l. Let the conditions (A. 4-) and {A. 5) hold true. Then, the system 
(A.l) must have a transition from (u, A) = (0, Aq), and there is a neighborhood 
U d X of u = such that the transition is one of the following three types: 

(1) Continuous Transition; there exists an open and dense set U\ C U 
such that for any (f Cz Ux, the solution ux{t,ip) of (A.l) satisfies 

lim \iinsvLp\\ux{t,Lp)\\x = 0- 



In particular, the attractor bifurcation of (A.l) at (0, Aq) is a continuous 
transition. 

(2) Jump Transition.- for any Ao < A < Ao + e with some e > 0, there is an 
open and dense set Ux d U such that for any ip G Ux, 

lim sup \\ux(t, ip)\\x > (5 > 0, 

t — >oo 

where S > is independent of X. This type of transition is also called the 
discontinuous transition. 

(3) Mixed Transition; for any Aq < A < Aq -I- e with some e > 0, U can be 
decomposed into two open sets and f/^ (Ul not necessarily connected): 

U=U^ + U^, n = 0, 

such that 

lim limsup ||M(t, iy9)||x = ytpEUi, 

^ ^-^0 t — *oo 

limsup \\u{t, ip)\\x > S > V(/J e C/^. 

t — >co 

An important aspect of the transition theory is to determine which of the three 
types of transitions given by Theorem |A.1| occurs in a specific problem. We refer 
the interested readers to [3J[T] for more discussions. Instead, here we consider the 
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transition of (A.l) from a simple critical eigenvalue. Let the eigenvalues f3j{X) of 
Lx satisfy 



(A.6) 




if A < Ao, 
if A = Ao, 
if A > Ao, 



Vj > 2, 



where /3i(A) is a real eigenvalue. 

Let ei(A) and el{X) be the eigenvectors of L\ and L*^ respectively corresponding 
to /?i(A) with 

^Aoei^O, L^^eJ = 0, <ei,ei>==l. 



Let ^{x, A) be the center manifold function of (A.l I near A = Aq. We assume that 

(A.7) < G{xei + <^{x, Ao), Ao), ei >= ax'^ + o{\x\''), 

where fc > 2 an integer and a 7^ a real number. 
We have the following transition theorems. 



Figure A.l. Topological structure of the mixing transition of 
(A.l I when fc=even and a ^ 0: (a) A < Ao; (b) A = Aq; (c) 
A > Ao. Here Ui is the unstable domain, and U2 the stable do- 



main. 



Theorem A. 2. Under the conditions {A.6) and {A.7), if k= even and a 7^ 0, then 
we have the following assertions: 



(1) {A.l) has a mixed transition from (0,Ao). More precisely, there exists a 
neighborhood U <Z X of u = such that U is separated into two disjoint 
open sets Ui and U2 by the stable manifold T\ of u = satisfying the 
following properties: 

(a) U = U^ + U^+Tx, 

(b) the transition in Ui is jump, and 

(c) the transition in U2 is continuous. The local transition structure is as 
shown in Figure A.l. 

(2) {A.l) bifurcates in U^~to a unique singular point on \ > Ao, which is 



an attractor such that for any ip G U2 , 

lim \\u{t,p)-v^\\x^Q, 

t — >oo 



where u(t,ip) is the solution of (A.l) 



(3) (A.l) bifurcates on \ < Xq to a unique saddle point v with the Morse index 
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(4) The bifurcated singular point can be expressed as 

= -(;3i(A)/a)i/('=-i)ei + oi\P,\'/^'-'^). 
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